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I. INTRODUCTION 


The improved characteristics of a modern missile make the interception process 
difficult: the reduced Radar Cross Section (RCS) makes it hard to detect by using 
a radar, while late emission leaves little time for intercepting and initiating any 
necessary defense and Electronic Counter Measures (ECM) activities. 

Conventional interception methods give only direction, and rough estimated 
range, and require a relatively long period of time. 

This work suggests a unique, non-conventional method for passive range and 
height estimation of a cruising missile, as well as any other microwave transmitter. 
The method is based on over sea multipath propagation, and uses five receiving 
antennas in a ladder configuration. For each receiving antenna there is a particu- 
lar multipath behavior, different than the other antennas. Comparing theoretical 
ratios of antenna pairs with actual received power ratios leads to an estimated 
location (range and height) of the transmitter. 

The motivation for the need of such a method is given in Chapter II. Chapter 
III describes the theory behind multipath propagation of electromagnetic waves, 
over different kinds of sea surface. It also explains in detail the model for rough 
sea which we used in this work. 

Chapter IV gives the basic theory of how the multipath propagation can be 
used for range estimation, and the conceptual estimation sequence. The detailed 
process, the assumptions being used during the research, and simulation results 


are given in Chapter V. 


Additional parameters that were not taken into account in the process, as well 
as those only mentioned, are considered in Chapter VI, along with their influence 
on the experiment. The results are analyzed. 

Chapter VII explains the computer program, and shows the concept behind 
the software process. 

Chapter VIII summarizes and concludes the whole work, and gives some sug- 
gestions for future steps. 

The rough sea model which we used in this work is the “Ament Model.” Other 
familiar models are described in Appendix A. Listing of the computer program used 


is given in Appendix B. 
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II. NATURE OF THE PROBLEM 


A. GENERAL 

In the modern battlefield, where the electromagnetic complexity and density, 
as well as a missile’s speed are very high, the probability of intercept is low, so 
detecting and intercepting missions became hard to achieve. The policy is to start 
the transmission as close to the target as possible, and to transmit for a minimum 
amount of time in order to be less “visible.” 

On the other hand, the radar cross section (RCS) of ships and especially 
missiles is being reduced, making it difficult for a radar to detect them. 

Therefore, we see that intercepting a missile, especially an ultrasonic one, is a 
tough mission and might be achieved only during the last seconds, which might be 
too late. 

By using a conventional interceptor, even if it is a fast one (monopulse re- 
ceiver), only a direction can be found. Using the ship’s radar for range finding can 
be worthless, because of the small RCS of the missile, and at the same time will 
help the other side to detect the source. 

For these reasons, we prefer to find a missile (as well as other ships’) location 
(direction and range) by using high speed passive equipment. 

Two common methods for solving this problem are being used, and are listed 
below: 

1. Triangulation 

This method requires two independent interceptors on two different plat- 


forms sharing the information of the target direction. The location of the target 


can be found on a map, as shown in Figure 2.1. Since each receiving antenna has 
its beamwidth, the position estimate is actually within a rhombus, and may not be 
as accurate as desired. The direction errors of the interceptors also add a particular 
level of inaccuracy to the estimation process. 

Other disadvantages of this method are the need for two ships to be in- 
volved, with all the logistics and synchronization, and obviously the time frame 


needed for this process, which make it a poor choice for missile defense. 


Real 
Target Location 


Area of 
Target 
Location 


> Estimated 






А 


аа c c 
Interceptor # 1 Interceptor # 2 


Figure 2.1: Triangulation Geometry. 


2. Auto-triangulation 


In order to estimate the location by using a single interceptor (single ship), 


we can make the triangulation by making sequential measurements, using the fact 


that the interceptor itself is moving and changing its place. Obviously, this method 
takes a long time for processing, and therefore is good only for stable targets (shore 
radar, for example). 

Actually, we cannot find the range of a moving target, especially a cruising 
missile, by using conventional methods, and so we need to have a special way for 
estimating the transmitter location. 

The method suggested in this work uses the multipath propagation of a 
transmitted electromagnetic signal in order to estimate the transmitter range from 
the receiver, while the direction is still being found using any common method such 
as rotating antenna, Instantaneous Direction Finder, etc. 

The theory of a general multipath propagation over smooth and rough sea 


is described in the next chapter. 


ПІ. THEORY OF MULTIPATH PROPAGATION 
OF ELECTROMAGNETIC WAVES 


In a usual radar case, the multipath propagation generally causes tracking 
difficulties, especially over a smooth sea, which causes specular reflections. 
In this chapter we will review the multipath propagation theory, divided into 


smooth and rough terrains. 


A. GENERAL 

A transmitted signal from a radar antenna spreads as it is leaving the antenna 
surface. A direct path exists between the transmitter (radar) antenna and the 
target (receiver antenna in our case), as well as an indirect path along which the 
signal transmitted from the transmitter antenna hits the sea surface and is reflected 
to the receiver antenna. 

The signals arriving from these two paths create interference in the receiver, 
and the relative phase difference dictates if it is a creative or destructive inter- 
ference. The path difference depends on the direct range R, transmitter antenna 
height À,, and receiver antenna height h,. Other factors have some influence as 
well, and will be discussed later in the next sections. Figure 3.1 shows the geometry 
of the multipath propagation over a smooth sea, assuming a flat earth. 

Since the sea state has significant contribution to the interference, a basic 
question that can be asked is how to distinguish between a smooth and rough 
surface (sea in our case). The answer for this question is not a straight-forward 


one, and, as a matter of fact, there is neither a single answer nor a defined border 
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Figure 3.1: Simplified Multipath Geometry Over Smooth Terrain and 
Flat Earth Assumption [From Ref. 1]. 


between the two cases. A method for having a quantitative solution for this problem 
was suggested by Lord Rayleigh, and is described herein [after Ref. 2]. 
Before starting the detailed analysis, we will give some definitions which are 


relevant to our discussion, referring to Figure 3.2: 
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Figure 3.2: Барај У of Irregular Surface and Path Difference [From 
Ref. 3}. 


e Wave Height Ah: The vertical distance between the wave’s trough and its 
crest. 


e Wavelength l: The horizontal distance between two adjacent crests (measured 
in the direction of wave propagation). 


e Wave Period T: Time interval between two adjacent crests passing the same 
fixed point relative to the ground. 


Consider an irregular surface with irregularities of height Ah. If two parallel 
rays of the same wavelength A hit the surface with a grazing angle v, as illustrated 


in Figure 3.2, the path difference between them will be Ar, where: 
Zx eas PN mn (3.1) 
Translating this into phase difference, we obtain: 


КОЕ БАТЕ Ar, 


and substituting (3.1) we have: 


Am AR 


Ag = \ 





sin y. (3.2) 


The phase difference Ag can vary between 0 and т, where a totally smooth 
surface will occur when Ag = 0. Zero phase can be achieved when there are no 
irregularities, i.e., Ah = 0, or when the grazing angle is zero. From the above, we 


can assume a smooth surface when: 


= +0 or № 0. (3.3) 


When the phase difference between the direct and reflected signal components 
is equal to z, there will be a destructive interference, and no received energy at all 


if there is a perfect reflection (|p| = 1). 


In order to choose a point where the surface becomes rough, we can arbitrarily 
choose any value between Aó = 0 and Aó = r. The most popular value is half the 
way (A@ = 7/2), as was chosen by Rayleigh, and hence, the “Rayleigh criterion” 


for a smooth surface is: 


À 
8sin ` 





Ah < (3.4) 


Other values have been chosen as well [Ref. 2,3,4,5,6]. 

From observing Equation (3.2), we can learn that the phase behavior is a 
periodic one, hence the interference will have maxima and minima in sequence. 

A typical received signal in the presence of multipath propagation over a 
smooth sea is shown in Figure 3.3. The amplitude decay is due to the spheri- 


cal spreading of the received signal (1/ R? . 
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Figure 3.3: Multipath Propagation. 


B. SMOOTH SEA 

Looking at the multipath propagation more carefully, we should notice that 
there are additional phase and amplitude changes in the reflected wave due to the 
reflection itself, as well as amplitude difference between the reflected wave and the 
direct one, since the propagation path lengths are not the same. 

The amplitude difference due to the difference in path lengths is very small, 
and hence can be neglected, since the paths differ only by a short distance compared 
to the total path length. [Ref. 7] 

The phase and the amplitude changes due to the reflection itself depend on 
the grazing angle, frequency, and polarization. 

Figures 3.4 and 3.5 illustrate the magnitude and phase (respectively) for verti- 
cal polarization, while Figure 3.6 shows the magnitude of the reflection coefficient 
for horizontal polarization. 

Note the small changes in the value of the reflection coefficient for horizontal 
polarization, which means that amplitude changes due to reflection can be neglected 
for small grazing angles (assuming reflection coefficient magnitude — 1 ). The phase 
difference in horizontal polarization can be considered as 180?, since for grazing 
angles of 0? to 90? it changes only up to 4°. [Ref. 7] 

Assuming a flat earth and a smooth terrain (reflection coefficient = —1), the 
ratio between the received power in the receiving antenna in the presence of mul- 
tipath, compared to a free space (for small grazing angles) is 7^ [Ref. 3,9], where: 


2лћ, (А, + В.) 


n? = 4sin? АСР 


(ф ¬+ 0), (3.5) 


using the same notation as in Figure 3.1. The shape of n* is shown in Figure 3.3 


for a typical case. 
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Figure 3.4: Reflection Coefficient in Vertical Polarization 
[From Ref. 7]. 


In order to be more precise, we need to take into account the fact that the 
earth is not flat. This fact will cause some changes in the geometry, and hence will 
influence the signal amplitude. 


Actually, the reflection coefficient depends on three main parameters: 


1) Plane smooth surface reflection coefficient p. 

2) Earth curvature, D,which reduces the reflection (in some cases, depending on 
the geometry). 

3) Surface roughness Ry. 


The actual reflection coefficient, hence, is [Ref. 6]: 


LEM (3.6) 


11 


“REET TTT we 
~ АМУ У О 
к МАҸ Г 22 00 
"АМАА ES 
А 

«ANNI 
„Т МУЫ 
m BEREIT. 


0 T" 2. = 4 50776 ои 
GRAZING ANGLE, DEGREES 


PHASE NEM DEGREES 


Figure 3.5: Phase Difference in Vertical Polarization [From Ref. 7]. 


The value of R, is between 0 and 1, depending on the terrain, and it will be 
discussed in detail in the next section. 

The earth curvature factor D (Divergence factor) has also values between 0 
and 1, where in the flat earth case it will be unity. 

Practically, in the scope of this work, we can use the flat earth assumption, 
since the geometry of the problem (transmitter height less than about 100 meter, 
receiver antenna height at about the same and both heights are much smaller than 
the distance in between) leads to small grazing angle and allows the assumption of 
parallel direct and reflected rays [Ref. 7]. 

Even though the transmitter 1s considered a point source, the reflections from 
the sea come from a particular zone, rather than in a single point. We shall look 


for all the points (loci), from which the reflected signal will reach the receiving 
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Figure 3.6: Reflection Coefficient in Horizontal Polarization 
[From Ref. 8]. 
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Figure 3.7: First Fresnel Zone [From Ref. 2]. 


antenna with the same phase shift, relative to the direct radiation, as illustrated 


in Figure 3.7. The loci obey Equation (3.7). 
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R, + R, — Ri = 6, (3.7) 


or after arranging: 


R. + R> = Ба + ó = Constant. (3.8) 


Equation (3.8) describes ellipsoids, so that the reflection zones are actually 
ellipsoids which form the Fresnel zones (Fresnel ellipsoids). Near the wavelength 
that has the most interest to us, 3 cm, (or frequencies around 10 Ghz), the area 
of the first Fresnel zone is about 2000 square meters [Ref. 2,3,6]. It is common to 


take the whole first Fresnel zone as the reflection area. 


C. ROUGH SEA 

The smooth sea case is generally rare, while some kind of roughness is more 
common, producing complications to the problem by adding more parameters, 
some of them arbitrary. Due to these facts, the rough sea case is difficult to 
analyze and predict. 

A more accurate solution leads to more complications and difficulties, so in 
practice we try to simplify and make assumptions in order to get a solvable case. 
Still we need to be very careful so that the result will approximate reality even 
after using these assumptions. 

In general, we can say that the first problem is to find a correct compromise 
between excess rigor which may yield no practical result, and too many simplifica- 
tions which may lead to incorrect results and hence wrong conclusions [Ref. 3]. 

In practice, the common general assumptions which take place in the rough 
sea case are: 


e Scattering elements are much larger than the propagated wavelength, 
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e Only the far field is calculated, 
e Diffraction and reflections between waves are neglected, 


A particular model for the roughness shape is assumed (depending on the 
case), and average wave height is assumed, and, 


e Constant sea state (roughness) is assumed for the whole area of interest. 
1. Reflection types 
At this point we will distinguish between two kinds of reflections: 
a. Specular reflection 

Similar to the smooth surface (as in a mirror case), specular reflec- 
tion can be characterized as obeying the classical optics rules, as well as having 
a coherent phase. The direct source for this kind of reflection is the first Fresnel 
zone (Fresnel ellipse) radiated points. The phases of the different points (from the 
receiver point of view) are approximately equal, hence the reflection is coherent, 
as mentioned previously. The specular reflection fluctuations are relatively small 
in amplitude. 

b. Diffuse reflection (scattering) 

Caused due to irregularities in the surface, scattering cannot be pre- 
dicted from the phase and directivity points of view. This scattering will occur 
over larger areas than the specular reflection (first Fresnel zone) and is incoherent 
in phase. The diffuse scattering has large amplitude fluctuations. We can also say 
that the smooth sea analysis is deterministic, while the rough sea analysis (both 
specular and diffuse reflections) is stochastic in nature [Ref. 5]. 

In general, there will be both reflection components, since a surface is 


neither totally smooth, nor extremely rough. The two components are illustrated 


in Figure 3.8 {Ref. 2,7,10]. 
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The specular and diffuse reflection coefficient components will be de- 
noted as p, and pa, respectively. The effective reflection coefficient will be a mul- 
tiplication of the smooth terrain coefficient pọ, and the relevant specular or diffuse 


coefficient. Hence, for the specular component, p,, the reflection coefficient will be: 


P = Po ` Ps, (3.9) 


and obviously will be smaller than that of a smooth sea po (since each coefficient 
is less than unity). 


The diffuse reflection will have a coefficient of 


р = ро · ра, (3.10) 


again, smaller than po. 
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Figure 3.8: Illustration of Diffuse and Specular Reflections 
[From Ref. 10]. 
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The diffuse reflection coefficient, due to its random nature, is found 
to have a Rayleigh distribution probability function for the amplitude, as given in 
Equation (3.11), while the phase has a uniform distribution between 0 and 27. 


The diffuse reflection coefficient amplitude has the form: 
pa = V2-T - exp(joa), (3.11) 
where: 


Og is the RMS wave slope (unitless), as shown in Figure 3.2, 


T = 0.77 - [1 — exp(—476)] - exp(—4.73 - 6), and, ЗЕ 
Oh ° SINY 


= 3.13 
б = 2000, (3.13) 


where øp is the standard deviation of the wave height, related to the peak to trough 
height. [Ref. 5] 

Assuming a smooth terrain, there will be a center point of reflection. 
The diffused reflection comes from the whole area surrounding this point, the glis- 
tening surface [Ref. 2,11], which depends on the observation angle from the trans- 
mitter to the receiver (E; in Figure 3.1), and on the RMS wave slope (o4). 

A typical RMS slope, og, is about 0.2 radians, and it depends on the 
wind velocity. The amplitude of the reflection coefficient pg has an extreme of 0.4, 
which rarely occurs, while a typical value is 0.35. [Ref. 11] 

We should notice that for a flat earth assumption, the grazing angle 
Vv 15 the same as the observation angle E, from the transmitter to the receiver 
(see Figure 3.1). The idea of the different reflection areas due to observation angle 
E,, and sea state (roughness) is demonstrated in Figure 3.9. In this Figure, 0, 


represents the rough surface RMS slope. 


l; 


The specular reflection coefficient has been analyzed by many, and 
there are quite a few models for this component. Each model uses different as- 
sumptions for the computations, and therefore there are (in some cases) different 
results. The conclusion from this is that the results are highly dependent on the 


roughness type which is assumed. 
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Figure 3.9: Reflection Areas in Different Roughness and Observation 
Angle, E.: (a) Specular reflection, p, ~ 1,p = po; (b) Slight 


roughness, (c4/À) - sin E, ~ 0.06, (c) Rough Surface [From 
Ref. 11]. 


In order to be as close to the experimental data with a generic formula- 


tion as possible, we can, in general, define the reflection coefficient, p, as [Ref. 12]: 


оо 


p= | р(Ф) соѕ Ф 40. (3.14) 


== © 


18 


where: 


ф 18 the phase of each particular element of the reflected wave in relation to a 
median plane surface, and, 


p(¢) is the probability distribution of all the phases ¢ along the path. 
2. Profile Types 
Using Equation (3.14), we find that the common profile models for rough 
terrain are (Ref. 2,12]: 
a. Sawtooth profile 
The phase ¢ varies uniformly between maximum values of +m, and 


the probability distribution therefore is: 





1 
m 3.15 
м9) = 51. (3.15) 
while the reflection coefficient is: 
DE а Pm = since om: (3.16) 


The idea here is that the wave height is uniform, but wave slopes are arbitrary. 
b. Sinewave profile 
The phase ¢ behaves like a sine wave with maximum amplitude of 


+ф,һ. In this case, the probability distribution is: 


1 
27ó cos u` 


р(Ф) = Gam 


where v is an arbitrary parameter, and: 


p Jo(óm ), 


where Jo is the Bessel function of zero order. 
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c. Random profile 
The phase ¢ has normal (gaussian) probability distribution with stan- 


dard deviation of ø. Here, 


_ exp — = Й 


= 3.18 
р(Ф) уру б ( ) 

and the reflection coefficient in this case 1s 
p = exp —(0* /2). (3.19) 


The latter profile is similar, in a sense, to the Ament model, which will be discussed 


later in this chapter. The three profiles are drawn in Figure 3.10 [Ref. 12]. 
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Figure 3.10: Reflection Coefficients of Common Models for Irregular 
Terrain [From Ref. 12]. 


A further discussion on other models is given in Appendix A. 

The most general and the most popular model is the one by W. S. 
Ament.* This model assumes a normal (gaussian) distribution, and it is popular 
since it does not include any discontinuities as other distributions do, but still 
meets the actual (experimental) data. For these reasons it seems to be the most 
realistic model (Ref. 2]. 

The Ament model is discussed and analyzed in many books and arti- 
cles, and from now on we will concentrate especially on it [Ref. 2, 3, 5, 6, 7, 9, 11, 


13, 14, 15]. 


D. AMENT MODEL FOR IRREGULAR TERRAIN 

The Ament model assumes that the wave height is normally (gaussian) dis- 
tributed, and it does not take into account the shadowing effect, where parts of the 
surface which are hidden by higher parts are not being illuminated by the incident 
wave. The shadowing effect takes place at very small grazing angles, while for 
higher angles it has no effect. The geometry is shown in Figure 3.11. [Ref. 13] 

Figure 3.12 shows the influence of the shadowing effect on the wave height 
standard deviation, related to the grazing angle [Ref. 13]. 


The specular component of the reflection coefficient in this model is: 


)3.20( ا 


3. == —2 
ps — exp ( n 


where: 
anp is the RMS deviation of the waves height, 


i 1s the grazing angle (refer to Figure 3.1), and, 


* W. S. Ament, “Toward a Theory of Reflection by a Rough Surface,” 
Proc. I.R.E., vol. 41, p. 142, January 1953. 
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Figure 3.11: Shadowing Effect Geometry [From Ref. 13]. 
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Figure 3.12: Shadowing Effect on Wave Height Standard Deviation 
[From Ref. 13]. 


A is the incident wavelength. 

Figure 3.13 shows the specular reflection coefficient component according to 
Ament model [Ref. 7]. 

The parameter øp is proportional to the sea state, and is given with other 


relevant parameters in Table 3.1 [after Ref. 2,11]. 
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Figure 3.13: Rough Surface Reflection Coefficient According to Ament 
Model. 


TABLE 3.1: WAVE HEIGHT PARAMETERS AND CRITICAL 
ANGLES FOR SPECULAR REFLECTION, VERSUS SEA STATE. 


sea state | description wave RMS height | critical angle (°), Ymax 


number of sea height (m) aon (m) 3 GHz 10 GHz 
0 calm 0 0 
1 smooth 0-0.3 0-0.065 » б > 1.8 
2 slight 0.3-1 0.065-0.21 1.8-6 0.5-1.8 
3 moderate 1-1.5 0.21-0.32 1.2-1.8 0.3-0.5 
4 rough 1.5-2.5 0.32-0.54 0.7-1.2 0.2-0.3 
5 very rough 2.5-4 0.54-0.86 0.4-0.7 0.12-0.2 
6 high 4—6 0:86-1.3 0.3-0.4 0.09—0.12 
(| very high > 6 О < 0.3 < 0.09 
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Figure 3.14: Scattering Regions Versus Roughness [From Ref. 11]. 
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Using the Ament model, we can distinguish more carefully between rough and 
smooth sea, and check the conditions for diffuse and specular reflection. Recall 
(3.20) for convenience: 


Ps = exp ( 


Aron sin Yh 
À 


We can consider the terrain as a rough one if p, < 0.5, which is equivalent to: 


Т sin Ymax > 0.0665. (3.21) 


For a particular wavelength at a given sea state, there is a critical grazing 
angle Umax, calculated from (3.21), above for which the sea surface is taken as 
rough. This angle is shown in Table 3.1. 

The specular reflection coefficient component is significant until the grazing 
angle exceeds twice the critical angle. Beyond that point, the diffuse component 
becomes the significant one. The different regions of roughness type are illustrated 


in Figure 3.14. [Ref. 2,11,14] 


Another way of looking at the specular (coherent) and diffuse (incoherent) 
scattering regions is glven in Figure 3.15. Here we can see the boundaries between 
the scattering regions for different grazing angles. The graph in Figure 3.15 is given 
for 2 GHz. For higher frequencies, the graph will become more compressed, by the 


reciprocal of the frequency. [Ref. 15] 
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Figure 3.15: Scattering Regions Envelopes at 2 GHz [From Ref. 15]. 


A similar graph for a constant wave height of 0.92 meter and different frequen- 
cles 1s given in Figure 3.16. Here again, for higher wave heights, the envelopes will 


be compressed by the same factor by which the height is increased. [Ref. 15] 
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Figure 3.16: Scattering Regions Envelopes at Peak to Through Wave 
Height of 0.92 meter [From Ref. 15]. 


IV. THEORY OF RANGE ESTIMATION 
USING MULTIPATH EFFECT 


A. THEORETICAL SOLUTION 

The multipath effect attenuation, as a function of range over a smooth sea, 1s 
well known, if both the transmitter and receiver height, as well as the frequency of 
the signal are given. This was described in the previous chapter. 

For a rough sea, the overall shape remains similar, except that the curve has 
more fluctuations. Still, by using a particular model, such as the Ament model in 
our case, we can deal with the multipath propagation as a known factor. 

Figures 4.1 and 4.2 show typical multipath propagation over a calm sea (state 
0) and over a very high sea (state 7), respectively. We can notice that in the rough 
sea case, the attenuation coefficient decreases arbitrarily in the short ranges, and 
has an arbitrary ripple over the whole range. 

From the above, the multipath attenuation 7^ can be written as a function of 
the range between the transmitter and the receiver RÊ and of the transmitter height 
hy, for a given receiver height hr: 


0^ — f(R, hi) | (4.1) 


h,.-const. 


The received signal 1s attenuated relative to free space. It 1s a combination of 
both the direct and indirect components. Hence, the received signal can be written 
as: 


Se == Sa : ^, (4.2) 


where: 

S, 1s the received signal power, 

Sq is the direct component of the signal power, and, 

л” is the multipath power attenuation factor, given by Equation (3.5). 

Since, as mentioned before, we cannot specify the direct component Sg, we 

need to eliminate it in order to solve the problem. A simple way for doing that is 
by using a pair of receiving antennas, one above the other, so that the range R will 
remain the same, while their heights have no practical influence on the direct signal 
amplitude. We should also remember that the transmitter height, h,, is constant 
for both receiving antennas. 


If we use the ratio of the output powers from the two antennas, 


(Sa = 54 · БЕН ni (4.3) 


(5:2 Sa'n n2 

we will be dependent only on the geometry of the receiver and the transmitter, 
without caring about the direct signals’ amplitude. 

Referring to Equation (4.1), we notice that we still have two unknowns, R and 
h,. Using the fact that we have two receiving antennas, while the range R and the 
transmitter height h, are the same, we need to find the correct ratio 72/792 which 
will give the right combination of range and height (R and h,). Figure 4.3 shows 
the geometry of the problem. A more detailed explanation of the sequence is given 
below. 

Let us assume, for simplicity, that the multipath propagation curve is a linear 
function (we shall return to the actual case later on), hence it behaves like an 


isomorphic function. 
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Figure 4.1: Multipath Propagation Over a Calm Sea. 
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Figure 4.2: Multipath Propagation Over a Very High Sea. 
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Figure 4.3: Range Estimation Basic Geometry. 


The power ratio of the outputs of the two receiving antennas can be written 


as: 








(Si ті Қа», ЕЗ | (4.4) 


(S)s п Ка “(Ве 

so it actually depends on the receiving antenna heights, which are known. 
Since we assumed that the problem has an isomorphic characteristic, a given 

ratio has only a single combination of range and transmitter height which fits this 

ratio. 


Keeping in mind those conclusions, the range estimation process will require: 


1) receiving a signal and analyzing its frequency and direction. 


2) using Equation (3.5), creating two (off line) tables (matrices) at that particular 
frequency, for which the coordinates are range and transmitter height, while 
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the contents are the attenuation coefficients for each R and h; combination, 
while the parameter is the receiver height (hence there are two tables). 


3) calculating an (off line) look-up table for the same frequency, ranges and trans- 
mitter heights, which has the values of the ratio between the two tables created 
in step 2. 


4) calculating the ratio between the two received signals, which is similiar (see 
Equation 4.3) to the ratio of the attenuation coefficients. 


Now, all we have to do is to search in the look up table from step 3, and find 
the range and the transmitter height, so the last step will be: 


5) searching for the right combination of range and height. 


B. PRACTICAL SOLUTION (OVER SMOOTH SEA) 

Practically, the multipath propagation leads to a multi-valued curve, hence 
the problem is not isomorphic any more, but has an ambiguous solution. Since 
the function is not linear any more, and involved with two unknown variables, 
we cannot use any analytical solution (like we could have used in the theoretical 
case), but only a numerical solution can take place. The idea is similar to the one 
described in the theoretical case with a pair of antennas. 

We saw before that two receiving antennas give a single ratio, hence we cannot 
solve the ambiguity problem (because of the periodic behavior of the multipath 
effect ). 

Adding a third antenna will give a second relation, so we can now use the 
process described previously for the two pairs, with a few candidate combinations 
(А and h,), and by comparing their locations between the two lists of candidates, 
choosing that one which appears in both of the lists; all the other locations are 


false ones. 


If the solution still is not a single one, we can add another antenna(s), until 
the comparison between all of the ratios leads to a single location. An important 
condition is that all of the receiving antennas be installed one above the other. 

A trade-off between the number of antennas and accuracy should be done: 
on one hand we want minimum ambiguous points (different points with the same 
ratios group) which dictates many antennas, but on the other hand we should use 
a realistic number of antennas. 

The influence of the different parameters on the multipath behavior can be 
seen in Figures 4.4 to 4.6. 

From observing these figures, we can conclude that: 


1) as the frequency increases, the curves become more complex (signal strength 
varies more rapidly with range), 


2) complexity increases with transmitter height, and, 


3) for increased receiving antenna height, the curves become more complex. This 
fact suggests keeping the antennas at lower height, but on the other hand, we 
should separate them enough to reduce the ambiguity. The receiving antennas 
in Figure 4.6 are 0.2 meters apart, so the curves are very similar to each other. 
In Figure 4.7, the antennas are separated by 3 meters, and the curves are well 
distinguished. 


Note that for short range we have very high complexity, which will be treated later 


on. 


C. PRACTICAL SOLUTION (OVER ROUGH SEA) 

In the case of rough sea, the main idea remains the same, although the signal 
exhibits fading as a function of the sea state. As the receiving antenna is raised, 
the influence of the roughness becomes more important. Figure 4.8 shows the 


multipath propagation curves over sea state 7, for different antenna heights. It 
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is obvious from these graphs that the higher the antenna is, the more rapid the 
fading. 

The look-up tables should take in account also the sea roughness parameter, 
where in our case we choose the Ament model, as described in Chapter III. 

After considering the facts mentioned above, the estimation process remains 


the same as described in the previous sections. 
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V. EXPERIMENTAL PROCEDURE AND RESULTS 


Now that we have presented the theory of how to use the multipath effect 


for estimating range, we can discuss how we applied that for solving the problem 


in our research. The assumptions, practical difficulties, and then some simulation 


results will be presented. 


A. 


ASSUMPTIONS 


Since the actual problem is a very complicated one, due to many surface 


parameters, where some of them are unknowns and/or very difficult to predict, we 


made some logical assumptions in order to obtain a solution. 


Polarization: the seeker was assumed to have horizontal polarization, hence 
the reflection coefficient was taken as unity in amplitude, with a 180° phase 
shift. The differences for vertical polarization will be discussed in Chapter VI. 


Sea state: the whole area of interest, 1.e., the whole range between the trans- 
mitter and the receiver was assumed to have the same sea state condition. 


Earth curvature: a plane earth was assumed. A spherical earth case is dis- 
cussed in Chapter VI. 


Atmosphere: a homogeneous atmosphere was assumed. 


Movement: the suggested method considers relative movement between the 
transmitter and the receiver. For homing missile, this is obvious. In other 
cases, such as estimating range between own ship and a stable transmitter 
(shore based radar, for example), we shall assume that the ship is moving 
relative to the radar. If there is no relative movement at all, the result can be 
incorrect. 


Geometry: It was assumed that the antenna separation, from the point of 
reflection, was much larger than the signal wavelength [Ref. 2]. Since we 
are dealing with microwave signals (from 2 GHz at least), this assumption 1s 
valid. (If this condition is not true, the reflection coefficients have some more 
complexity.) 
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e Range interval: almost any range interval can be assumed, but note that 
for the long ranges we should check the geometry and see if the plane earth 
assumption is still applicable. Since we concentrated in this work on homing 
missiles, the seeker transmitting range is the one of importance, and hence a 
maximum range of 20 Km (about 11 NM) is considered, while a minimum of 
1 Km (about half a NM) assumed. 

B. ESTIMATION PROCESS 
With the assumptions of the previous section in mind, let us repeat briefly the 
range estimrtion process. We shall look more closely at the details shortly. 
1. General Discussion: 
A general solution requires the following steps: 
a. Frequency Measurement: 

Since the whole process is frequency dependent, this will be the very 

first step made by any EW receiver. 
b. Off Line Tables: 

In this step we build the theoretical tables for all ranges and heights of 
interest (i.e., the possible locations of the transmitter), for each receiving antenna 
height. The inputs for creating those tables are frequency, receiving antenna height, 
and sea state, which dictates the attenuation factor to be placed in the table. 

c. Attenuation Coefficient Ratio Tables: 

After we have the theoretical receiving values (from the previous step), 
we calculate (for each range and height combination ) the ratios of the output signals 
for each receiving antenna pair. 

d. Amplitude Measurement: 
We must measure the received signal amplitude of each receiving an- 


tenna (we will discuss later the number of receiving antennas), and make a table of 


amplitude ratios. For example: suppose we have three receiving antennas, where 
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each one receives a signal with amplitude $1,52, and s3, respectively. The ratio 
table will have two elements: s2/s; and 53/31. 
e. Comparison: 

Now, that we have the measured signal ratio (from step d), and the 
theoretical one (step b), we need to compare - for each receiving antenna pair, the 
actual value with the off line (theoretical) value, and record range-height combina- 
tions for which they are equal. These range-height combinations represent possible 
iransmitter locations. 

Since, as we mentioned in Chapter IV, the problem has an ambiguous 
solution, there will be several possible locations, even though only one of them is 
the correct one. The next step is to determine the correct location. 

f. Validation: 

In order to choose the correct location from all the candidate locations 
determined in step e, we will compare all the locations determined for the first pair 
of receiving antennas, and the second pair of receiving antennas. There is, ideally, 
only a single location, which should appear common to all the pairs. 

Practically, because of the high ambiguity level (especially for the rel- 
atively short ranges), and also because of approximations in the implementation, 
there will still be more than a single location. In order to resolve this problem, 
we shall use more receiving antennas, which provides more ratios. This allows the 
ambiguity to be reduced to an acceptable minimum. 

g. Output: 
The matrix location is translated into an estimate of transmitter range 


and height. This is accomplished by using the coordinates of the look up table. 


2. Detailed Discussion 
The following section explains in detail every step listed before. 
a. Frequency Measurement: 

Since this requires an external system, we will not discuss it. Any EW 
receiver on board can be used for frequency measurement and direction finding. 
The frequency data is provided as an input to the range estimation process. 

b. Off Line Tables: 

This step uses the theoretical equations for calculating the multipath 
attenuation factor n?. For each receiving antenna we have a particular table, where 
the coordinates are transmitter height, h: (vertical axis), and range, R (horizontal 
axis). 

The height axis is divided into 36 “slices” with a 1 meter resolution, 
starting at a minimum height of 5 meters, and ending at a maximum of 40 meters. 
(These numbers are parameters, which can be changed easily. ) 

The range axis runs from 1 Km to 20.5 Km, and was originally divided 
into 130 equal steps of 50 meters each. Since the ambiguity for short ranges is 
very high, this partition caused some problems in finding the correct location for 
short ranges. This suggested the need for higher resolution at short range. The 
chosen partition is 130 steps of 50 meters each, up to the range of 7.5 Km, and 
then reduced resolution of 125 steps, 100 meter each, until the maximum range. 
(Again, all these numbers are parameters which can be changed.) 

The multipath attenuation factor is calculated by using the following 


formulae: 


The basic attenuation, n^, for a smooth sea case is [Ref. 3,9]: 


4 2 2rh,(h, ae Ше) 


n? = 4 sın ХБ (ф — 0), (5.1) 
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where: h, 1s the transmitting aU NO height, 
h, is the receiver antenna height, 
A is the signal wavelength, 
R is the range between transmitter and receiver, and, 
ф 18 the grazing angle. 
In the case of a rough sea, we have to multiply 72? Бу the reflection 
coefficient, p, which is different in various models. In our case we used, as described 


in Chapter III, the Ament model: 


Этоь ant) (5.2) 


p = exp —2 ( T 
where: op is the RMS deviation of the waves’ height, 
w is the grazing angle, and, 
A is the signal wavelength. 
Therefore, in a rough sea case, the multipath attenuation factor will 
be: 
Поезда = Таан Р. (5.3) 
The grazing angle w in equation (5.2) is a derivative of the geometry. 
It depends on transmitter height, receiver height and range. Figure 5.1 shows the 
geometry from which we calculate w. 
We shall now derive the grazing angle ? [According to Ref. 7]. From 


Figure 5.1 we see that 
Ra = dc ee (5.4) 


We define the range difference between the indirect path and the direct 


path as ó, where: 
ó = В, + Г» = Ra (929) 
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E Earth surface 


њу 
N 


Figure 5.1: Flat Earth Multipath Geometry [From Ref. 7]. 


Looking at the image of the transmit antenna, we see that 
Ву + В; = V(hi + hy)? R2. (5.6) 
Substituting Equations (5.4) and (5.6) into Equation (5.5), we obtain: 
(hı, + hr)? + R* — J(h, — h,)* 4 R*. (5.7) 


Let us assume that the range R is much greater than the antenna 
heights, and also R? > (hı + h,)2. In this case, we can approximate Equations 
(5.4) and (5.6) as: 


(hi +h) | 
Ra = МБ е В |1 + Се ў (5.8) 
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and, 


h, —h,)? 
Ry + Ro = V (h, + Az)? = ІК 1+ | ; (5.9) 


respectively. 
Substituting Equations (5.8) and (5.9) into Equation (5.5), we obtain, 


after arranging: 


0 (hı + hr)” — (hı — hr) = 2h,h,. 


1 
2R R 0 
Observing again the imaged antenna, the angle 0, is given by 
h, + h 
0, — sin! —— —— 5.11 
ш Ку + R2’ (5-11) 


and we notice that it equals (interchangeable angles) the grazing angle, 2. Using 
Equation (5.5), we finally find the grazing angle to be: 


--1 he + Rh, 


ф = sin Ras" 


(5.12) 


Since we know all the parameters from solving previous equations, we 
can now find vy. (Actually, we do not know the range R, but we should remember 
that in this step we are building tables, so the range is being used as a parameter.) 

The only variable still missing for solving Equation (5.2), (the reflection 
coefficient), is the wave height RMS deviation øp. This value depends on the sea 
state, and is given in tables in the literature. Table 3.1 shows this parameter (as 
well as some others), as a function of sea state. 

The propagation attenuation factor ( 2 in the one way case) also 
has its impact on the received signal. Since we deal with ratios, this factor cancels 


itself, and therefore had not been taken into account during the calculations. 
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с. Attenuation Coefficient Ratio Tables (17 /5;): 

Having created the different tables for each antenna, we now create a 
new set of tables, which contain the ratios of the received signal amplitudes for 
the antenna pairs. For three receiving antennas we get two tables, while for N 
receiving antennas there will be N — 1 tables. 

Figure 5.2 shows graphs of the ratio table data for sea state zero, 
at 9 GHz and a transmitter height of 20 meters. Note that for clarity only one 
“slice” is shown, referred to the particular transmitter height. The table is actually 
three-dimensional (range, transmitter height, and amplitude ratio). 

d. Amplitude Measurement: 

The first and major assumption is that all the receiving channels are 
calibrated. In this case, two similar signal amplitudes at the input to different 
channels will create the same amplitude at the outputs of those channels. 

Actually, there is no need for single channel amplitude measurements, 
since we deal with ratios. For this reason, we measure the individual signal ampli- 
tudes and immediately divide between pairs to form the desired ratios. 

e. Comparison: 

We now have in hand all the details needed to find the correct location 
of the transmitter, assuming of course, that the model is accurate. 

A comparison is made among the ratio levels from steps c and d. Every 
place that both levels are the same, is “marked down” by a “candidate location 


22 


flag.” Figure 5.3 shows the received signal ratios (constant level line) and the 
corresponding candidate location flags, for four antenna pairs. 
Again, we should remember that the real comparison is being made 


between many more points, since there is a third dimension (height). For simplicity, 


we show here only one of 36 height slices. 
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The dots in Figure 5.3 (upper portion) represent the real calculated 
values, while the connecting curves are made by extrapolation. Therefore, a can- 
didate location flag is set only when the measured signal ratio is within a specified 
distance from the calculated signal ratio. This distance is determined by the width 
of a comparison window which is a software parameter. The wider the window, the 
more candidate location flags we will set. On the other hand, we cannot allow too 
narrow a window, since there is a risk of loosing correct locations. It is obvious that 
as we use higher resolution in the lookup table, we can also narrow the comparison 
window. 

Figure 5.4 shows the candidate location flags for the same scenario as 
in Figure 5.3, but with a wider comparison window. Notice that more flags appear 
in this case. 

For a higher sea state, there are random fluctuations in signal level, 
and the comparison becomes harder. We must open the window wider, but at the 
same time more incorrect locations are flagged. A practical compromise is required. 

f. Validation: 

After we have all of the candidate location flags, we need to decide 
which location is correct. The method of doing this is by doing a cross location: 
comparing between all the flags, and checking which one repeats in all of the ratios. 

If after this process there is more than a single location (more than 
only one cross location), the measurement cannot give a correct location, so we 
discard it, and start a new session. 

Practically, using three receiving antennas (two ratios) gives no single 
location in most cases, and therefore we need to add more antennas in order to have 
more ratios. During the experiments we tried using different numbers of antennas, 


and the final sufficient result uses five receiving antennas. 
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Figure 5.5 shows the candidate location flags for each antenna pair and 
the crossing between them, giving a single valid location marked by a diamond. The 
sea state is zero, and the window width is 2.5%. Here, we can see not only a single 
“slice”, but the whole combination of range and height. From this it is obvious that 
a single pair of antennas gives no information at all (too many candidate location 
flags). 

For the case of three antennas, (see Figure 5.6), the ambiguity is re- 
duced significantly (only two flags overlapping), while adding a fourth antenna, 
reduces the location to a single, correct one. In this particular case, four receiving 
antennas give the correct location. In many cases, though, four antennas will still 
leave some order of ambiguity, but the fifth one will, in most cases, give a correct, 
single location. 

The ambiguity level reduction as a function of number of receiving 
antennas is well illustrated in Figure 5.6, for sea state zero. 

For a wider window there will be more flags, as shown in Figure 5.7, 
and here we can see that only the fifth antenna solves the ambiguity. 

When the sea state is higher, a wide window will lead to a large level 
of ambiguity. Figure 5.8 shows a 12% window in sea state 3, where there are four 
locations. Reducing the width to 5% will reduce the number of locations to two, as 
shown in Figure 5.9, while the ambiguity will totally disappear for a 2.57% window, 
as in Figure 5.10. 

g. Output: 

The previous step brought up a single location flag, which shows the 
correct location of the transmitter in the range-height table. The output step 
translates the table location into the real coordinates of range and height, and 


types it as an output result. 
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C. FINAL RESULTS 
As mentioned earlier in this chapter, the final resolution used here, is 130 steps 
of 50 meters each, from the minimum range up to 7.5 Km, and then 125 steps of 
100 meters each to the maximum range. The comparison window width is 2.5%. 
With these parameters, estimation results are given in Figures 5.11 and 5.12, 
for sea states 0 and 3, respectively. It should be mentioned that the estimated 
range sometimes deviates from the actual range by a few tens of meters, which 


cannot be seen in the graph, due to resolution. 
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Figure 5.3: Comparison Process with 2.5% Window. 
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Figure 5.3 (continued): Comparison Process with 2.5% Window. 
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Figure 5.4: Comparison Process with 12% Window. 
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Figure 5.4 (continued): Comparison Process with 12% Window. 
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Figure 5.6: Ambiguity Level vs. Number of Receiving Antennas; Sea 
State 0 and 2.596 Window. 
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Figure 5.7: Ambiguity Level vs. Number of Receiving Antennas; Sea 
State 0 and 12% Window. 
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Figure 5.8: Ambiguity Level vs. Number of Receiving Antennas; Sea 
State 3 and 12% Window. 
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CANDIDATE LOCATION FLAGS FOR 3 PERCENT WINDOW 


Figure 5.9: Ambiguity Level vs. 
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Figure 5.10: Ambiguity Level vs. Number of Receiving Antennas; Sea 
State 3 and 2.5% Window. 
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Figure 5.11: Estimated Results (Range); Sea State 0. 
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Figure 5.11 (continued): Estimated Results (Height); Sea State 0. 
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Figure 5.12: Estimated Results (Range); Sea State 3. 
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Figure 5.12 (continued): Estimated Results (Height); Sea State 3. 
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VI. ADDITIONAL PARAMETERS CONSIDERATION 


This chapter deals with some parameters that have influence on the calcula- 
tions and terms used for the estimation process. A few of them were mentioned 


but not discussed for all cases, and the others were not discussed previously. 


A. VERTICAL POLARIZATION 
So far we dealt with horizontal polarization, and all the calculations were done 
for this polarization. Vertical polarization complicates the problem a bit, since the 
reflection coefficient, p, changes significantly with the grazing angle, as shown in 
Figure 3.4. 
Let us check some representative cases in order to see how much the reflection 
coefficient really changes due to vertical polarization. 
Case 1: transmitter is located 10 Km apart from the receiver, at a height of 
10 m, and it transmits with a frequency of 10 GHz. The receiver antennas are at 
heights of 12 m and 20 m. So, our case has the parameters: 
i = 10 km, 
ni l0 шї, 
J = 10 Gl (AS 3 стир 
la tlm and, 
ha = 20 m. 
Using equations (5.4) to (5.12), we find that the grazing angles are equal to: 


фу = 0.126?, and y? = 0.172°, respectively for the receiving antenna heights. 
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Using the correct graph, (Figure 3.4), we find that the reflection coeficient due to 
the reflection, equals 0.98 for both grazing angles. 

Case 2: in the second case the transmitter is higher, and the range between 
the receiver to the transmitter is reduced: 

R = 4 km, 

h, = 40 m, 

DEC HZ = 3 cm), 
hr = 12 m, and, 

hm = 20 m. 

Using the same method as in the first case, we find here, that the reflection 
coefficient for the lower receiving antenna is 0.82, while for the higher antenna it 
is 0.80. 

Case 3: here we take an extreme case of very high transmitter, as well as 


receiver antennas: 


ic = 4 km, 

h, = 100 m, 

ЕОС На (А = 3 ст), 
ћа = 38 m, and, 

hy. = 50 m. 


In this case, the reflection coefficients are 0.58 for the 38 m receiving antenna, 
and 0.55 for the 50 m antenna. 

We can conclude that for low geometry (transmitter and receiver), and large 
ranges, the reflection coefficient can be assumed the same for all receiving antennas. 
For shorter ranges and larger heights, we should take into account the changes in 


the reflection coefficient for the different receiving antennas. 
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Practically, for the vertical polarization case, we need to modify the off line 
tables by the amount of the reflection coefficient changes. This can be done, since 
we know how to calculate it with the given parameters, and for the whole range 
and transmitter height field. The only parameter that has doubt is the kind of 
polarization. This verification must be done by another system, and provided as 


an input to the estimation process. 


B. RECEIVER HEIGHT 

As we saw in Chapter 4, the higher the receiver antenna is, the more complex 
the curve is. In our work we used a maximum receiving antenna height of 20 meters, 
since it is a reasonable height for a conventional missile boat. For a larger ship (such 
as a corvette or destroyer, for example), we shall need higher receiving antennas. 
The only influence, as mentioned before, is on the curves complexity. The solution 
for this is to use a greater resolution in the look up es to compensate for the 


higher ambiguity. 


C. SHIP TILTING 

For the ships of our interest (relatively small missile boats), the worst case 
rolling angle can be around 15° for cross wind cruising direction in sea state 3. 
For sea state 4, the extreme angle is about 20° for the same conditions. Since 
the tilt occurs for both sides, we can say that, to a first degree of approximation, 
the average tilt is 0°. Even though, let us check the worst case. The geometry of 
the problem is given in Figure 6.1. We will now check the tilt impact on the look 
up tables for two cases: the closest antennas (12 meters and 14 meters), and the 


extreme antennas (12 meters and 20 meters). 
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Figure 6.1: Tilted Ship Geometry. 

Referring to Figure 6.1, in our case, h; 2 12 m, hs = 14 m and hs = 20 m. 
The maximum tilt (roll) angle, “max, is taken as 15? for sea state 3, and 20? for 
sea state 4. 

From Figure 6.1 we see that the actual height of each receiving antenna, after 
being tilted, is: 

DES cosa G = 2..235. (6.1) 


For sea state 3 (Qmax = 15°), we obtain: 
осо = бт (where h, = 12 т), (6.2) 


and 


h; = 20 - cos 15° = 19.3m (where h} = 20 m). (6.3) 


The actual distance between the antennas is 7.7 meters, instead of a theoretical 
value of 8 m. Figure 6.2 shows the ratio 52/5? for receiving antennas located 


perpendicular to the sea level, in heights of 20 meter and 12 meter, respectively. 
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Figure 6.2: Ratio Between Attenuation Coefficients; Sea State 3, 
Receiving Antennas at 12 m and 20 m. 
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On the same figure we see also (n2)'/(n£)' for the tilted case. Figure 6.3 shows the 


ratio between those ratios: 
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Figure 6.3: Ratio Between Ratios of Attenuation Coefficients; Sea State 
3, Receiving Antennas at 12 m and 20 m. 


Now we will repeat the calculations for receiving antennas at heights of 14 


meters and 12 meters. Note that hj remains the same (11.6 m), while ^ becomes: 


E coll 1305 (where h> = 14 m). (6.4) 
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Figure 6.4: Ratio Between Ratios of Attenuation Coefficients; Sea State 
3, Receiving Antennas at 12 m and 14 m. 


The actual distance here is 1.9 meters, instead of 2 meters. Figure 6.4 shows 
the ratio between 72/5? (perpendicular antennas), and (7n2)'/(n?)! (tilted anten- 
nas). 


Repeating the same calculations for sea state 4 (@max = 20°), we obtain: 
h = сас ценне m: (where h, = 12 m), (6.5) 


and 


h = 20۰cos20° = 18.8 m (where hs = 20 m). (6.6) 


~] 
to 


The actual distance between the antennas is 7.5 meters, instead of a theoretical 
value of 8 meters. Figure 6.5 shows both ratios of 72/n? and (n3)*/(n})*, while 
Figure 6.6 shows the ratio between those two ratios. 


For the closest antennas we obtain: 
h; = 14 : cos 20° = 13.2 m (where h; = 14 m), (6.7) 


while h) remains the same, 11.3 m. The actual difference is, therefore, 1.9 meters, 
instead of 2 meters. Figure 6.7 shows the ratio between 52/5; (perpendicular 
antennas), and (n4)? /(n1)? (tilted antennas). 

To summarize these results, we see that the method is sensitive to ship tilting, 
since there is a difference between the tables for perpendicular and tilted config- 
urations. The conclusion of these results is that the receiving antennas must be 


stabilized. 


D. RECEIVER NOISE 

In the case of our work, the received signal to noise ratio was assumed to be 
relatively high. The ranges of interest are less than 20 Km, while the missile ERP 
is in the order of 100 to 120 dBm. In addition, the transmitted signal propagates 
only one way. Therefore, any received signal is much larger than the receiver noise, 


and the signal to noise ratio is high enough, so we can disregard the noise. 


E. ATMOSPHERIC ATTENUATION 

The atmospheric attenuation can be treated in two ways: modify the look up 
tables with the attenuation data, since all the relevant parameters are known, or, 
more practically, to disregard it. Since we deal with ratios, the attenuation will 


have an influence only if the path difference is significant. As we already discussed 
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Figure 6.5: Ratio Between Attenuation Coefficients; Sea State 4, 
Receiving Antennas at 12 m and 20 m. 
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Figure 6.6: Ratio Between Ratios of Attenuation Coefficients; Sea State 
4, Receiving Antennas at 12 m and 20 m. 


in Chapter III, the path difference is negligible compared to the total path length. 
Therefore, the attenuation is practically the same for the direct and indirect path, 


so it has no influence on the problem, hence can be neglected. 


F. OTHER MODELS 

Appendix A gives other models for reflection coefficients, instead of the Ament 
model, that we used in this work. To use another model, the only thing that should 
be taken care of is modifying the model equation in the computer program. This 


point will be discussed in detail in the next chapter. 
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Figure 6.7: Ratio Between Attenuation Coefficients; Sea State 4, 
Receiving Antennas at 12 m and 14 m. 


G. EARTH CURVATURE 

In our work, we assumed a plane earth. In general, some geometries need 
to take account of the earth’s curvature. Our particular geometry results in no 
influence from the earth curvature, as shown below. The extreme condition 15 for 
the lowest receiving antenna, and the lowest missile. 

The maximum range that still allows plane earth assumption is calculated by 


using equation 6.8 [Ref. 7]. The problem geometry is shown in Figure 6.8. 


Rmax = \/2а„Һ„-+Е h2 /2ach Thr, (6.8) 


| 


Figure 6.8: Plane Earth Assumption Geometry. 


where: 

a, is the effective earth radius (= 4/3 x В), 

R, is the earth radius (= 6370 Km), 

h, is the receiving antenna height, and, 

h, is the transmitter height. 

For all practical purposes, the term of h^ under the square root can be ne- 
glected. 

Our extreme case is when the missile cruises in hy = 5 m, and the lowest 
receiving antenna is located at height of hy = 12 m. For this geometry we obtain 


a maximum range of: 


Hea V2 5493-00): 12 V2 - 8493-10" - 5 = 23.5 Km. 


ТТ 


Since we deal with relative ranges of up to 20 Km, we can | 


curvature. 
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VII. COMPUTER PROGRAM 


The computer program in this work was written in APL. The listing is given in 
Appendix B and it includes all the necessary comments. This chapter presents an 
overview of the whole program, and then explains each function in greater detail. 

The program consists of a few functions which call each other in sequence, and 
a main management function. 

The sequence follows exactly after the one described in Chapter V, under 
“Detailed Discussion.” All the inputs that must be given by external means, are 
inserted here by the user, as software parameters. 

The program is started by typing RUN. The user is asked to enter sea state 
(0 to 7), and with this data, the RUN function calls the main function for receiving 
data, named DATAIN. This function automatically creates more needed parame- 
ters, in a random manner, between given limits. It is also possible for the user to 
dictate each one of these parameters, instead of creating them randomly by the 
program: 


е frequency — between 2 to 18 GHz, 


e ERP (Effective Radiated Power) - between 100 to 120 dBm (10^ to 10? 
Watt), 


e initial range for simulation (range where missile starts to radiate) — between 
10 to 20 Kilometers, and, 


e transmitter height — between 5 to 40 meters. 
Other relevant parameters are programmed off-line as software parameters, 


and can be changed easily at any time. Those are: 
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e window width, as discussed in Chapter V, 

¢ minimum transmitter height for the look up table, 

e number of steps for the transmitter height for the table, 
e resolution of transmitter height steps for the table, 


e minimum range for simulation (if the missile passed this range without being 
intercepted, the simulation is failed and stops), 


e range resolution for the look up table, both high resolution for short range, 
and regular resolution for longer range, 


e number of range steps for each range resolution, 
e resolution of simulation steps (range distance between each receiving), and, 


e locations (heights) of receiving antennas. 

The function then calculates the maximum range, and shows the major pa- 
rameters on the screen (a full session is demonstrated later in this chapter. ) 

With all the necessary data, DATAIN returns to the management function, 
RUN. In the next step, RUN sends all the relevant parameters to the main function 
STARTPRO, and activates it. 

STARTPRO is responsible for all the major calculations: 


е creating off line tables which include also the attenuation coefficient, 77, 
and the reflection coefficient, p, 


e receiving five signals in the receiving antennas, and divide their powers to 
determine the ratios, 


e comparing between the off line tables and the actual ratios, 
е marking the correct (single) location, and finally, 


e sending the location to the inanagement function. 


S0 


The signals for STARTPRO are created by the simulator SIM. This function 
is activated by STARTPRO, every time it needs a new set of receiving signals. 

SIM receives the relevant parameters from DATAIN, and creates five signals 
due to those parameters and the geometry. The signal is calculated using equa- 
tion 7.1. 


S;= ERP -1*-p- (7.1) 


m 
4r R?’ 
where: 

ERP is the effective radiated power of the transmitter, 

R is the relative range between the transmitter and the receiver, 

n^ is the attenuation coefficient, calculated by using Equation (3.5), and, 

p is the reflection coefficient, calculated by using Equation (3.20). 

After the management function RUN received the location of the transmitter 
from STARTPRO, it calls the translating function RHFIND in order to translate 
the coordinates of the table into range and height, and to type it. 

A conceptual flow chart is given in Figure 7.1, and a simulation session is 


shown in Figure 7.2. 
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Figure 7.1: Estimation Process; Conceptual Flow Chart. 


Figure 7.2: 


* RUN 
ENTER SEA STATE: 


0 - CALM ] - SMOOTH 2 - SLIGHT 
3 - MODERATE 4 - ROUGH 5 - VERY ROUGH 
6 - HIGH 7 - VERY HIGH 


*2 (for example) 


RANGE INTERVAL FOR SIMULATION: 
MINIMUM RANGE (IN METERS) 1000 
MAXIMUM RANGE(IN METERS) 20000 

SEA STATE: 

2 

FREQUENCY (IN GHZ): 

8.85 

INITIAL RANGE (IN KILOMETERS): 

14.5 

TRANSMITTER HEIGHT (IN METERS) - ACTUAL: 
36 

RECEIVER ANTENNAS HEIGHTS (IN METERS): 
(ea 45 1195886250 


TRANSMITTER RANGE (IN KILOMETERS): 
14.4 

TRANSMITTER HEIGHT (IN METERS): 

36 


UPDATED ACTUAL RANGE (IN KILOMETERS): 
14.4 


Simulation Session Example. (Lines Starting with Asterisk 
(*) are Entered by the User.) 
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VIII. SUMMARY AND CONCLUSIONS 


The work which had been done during this thesis research, together with the 
results obtained, show that the suggested method for passive range estimation 
using multipath effect can probably be used. The work developed the method 
theoretically, and proved its feasibility. 

The accuracy of the estimated range and height depends upon many param- 
eters, as well as receiving antennas configuration. A trade off should be done be- 
tween number of receiving antennas, their heights, separation and window width, 
as described in Chapter V. 

There is no analytical way for finding the optimum trade off, and only trial 
and error sessions can achieve a local optimum, as was done during this research. 

Another question of interest is the accuracy of the rough sea model (Ament 
model in our case). 

The answer to this question can only be obtained by sea trials. Those trials 
should answer the questions listed below: 


e Is the theoretical analysis of the additional parameters in Chapter VI cor- 
rect, or should we modify these parameters? 


e What is the estimation accuracy if an incorrect sea state is entered, and 
is there a preferred direction for entering sea state (is it better to enter a 
higher sea state in case of uncertainty, or lower state)? 


e How does the wind direction influence the accuracy, if at all, (upwind, down- 
wind, and crosswind)? 


e What is the practical feasibility of the suggested method, and with which 
limitations (such as sea state and wind direction)? 
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We should notice that the Ament model neglects the shadowing effect, al- 
though it deals with small grazing angles, where this effect does occur. It is ques- 
tionable how far the shadowing effect changes the reality from the model. 

A possible influence of the shadowing effect is the fact that the amplitude of 
the reflected wave will be lower, so that the signal minima will not be so sharp as 
is shown in Figure 4.2. A decrease in the amplitude of the specular reflection will 
also occur occur because of high sea state, and capillary waves. If sea trials show 
consistent errors, we can conclude that the Ament model is inadequate. In this 
case it will be necessary to develop another more accurate model. 

It can be predicted that for sea state 0, the results will be more accurate than 
in any other case, while the more roughness, the less accurate the results. 

In order to improve the estimation process, and to reach a lower rate of errors, 
we can modify the program. The modified program will consider also the past 
readings, and compare each present estimate with the previous ones. If the range 
becomes shorter by a reasonable increment, and the transmitter height remains 
almost the same, we can conclude that the estimate is correct. On the other hand, 
if the estimate does not meet this test, then another estimate should be made, and 


compared, and so on. 
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APPENDIX A: REFLECTION COEFFICIENT MODELS 


In Chapter III we discussed the most popular model for reflection coefficient 
of a rough sea, the Ament model. There are several other models, which are given 
below. It should be noted that there are additional models which we will not cover 


here. 


A. BULLINGTON MODEL 
This model distinguishes between different surface profiles [Ref. 3,12]: 
1. Sawtooth Profile 
The wave height is distributed uniformly, but does not necessarily have 


uniform slope. The phase values have equal probability between —¢, and +@m, 


where 
2ro, sinw 
Om ~ =y (A.1) 
and 
сь is the RMS deviation of the waves height, 
w is the grazing angle, and, 
A is the signal wavelength. 
The reflection coefficient for the sawtooth profile is given by: 
DI = =зшео (А.2) 


and is shown in Figure А.1. 
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Figure A.1: Bullington Model — Sawtooth Profile. 


2. Sinusoidal Profile 
The reflection coefficient in this case is given by Equation (A.3), and is 


shown in Figure À 2. 


p= Jo (у (А.3) 


where Jo is the Bessel function of the first kind and zero order. 


3. Triangular Profile: 


Here, the reflection coefficient is given by: 


_ f sin ó, 4 | 
p= ( ф, ) 2 (5.3) 
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Figure A.2: Bullington Model — Sinusoidal Profile. 
where the maximum phase is: 


лаћ ап 


Qe S REN (4.5) 


Figures A.3 and A.4 show the triangular profile reflection coefficient. 


4. Sine Square Profile 


The reflection coefficient is more complex in this model, and is given by 


Equation (A.6), and shown in Figure A.5. 


p = cos $1 - Јо(Ф,), (4.6) 


where ¢; is the same as given in Equation (A.5). 
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Figure A.3: Bullington Model — Triangular Profile. 
B. DAVIES MODEL 

The Davies model assumes that the surface has a random behavior, which 
yields a gaussian height distribution and correlation coefficient. This model neglects 
the shadowing effect, and is good for small to medium grazing angles (y < 40°). In 
order to describe the surface, we introduce an arbitrary parameter, m/o. Experi- 
ments above the sea surface show that a value of m/o = 50 to 60, gives sufficient 
correlation between theoretical and experimental results. [Ref. 3]. 


The reflection coefficient for this model is given by: 
т | т 2 2 = 
p = ç ` V seci : exp - 52) tan 1 (А.Т) 
Figure À.6 shows the reflection coefficient. The value m/2o is 27 (m/o = 54). 
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Figure A.4: Bullington Model — Triangular Profile (expanded). 


C. ISAKOVICH MODEL 

In this model we assume that the surface is irregular, isotropic in all directions, 
and has a normal height distribution. The parameter m/c describes the surface (as 
in the Davies model). It is the only parameter which dictates the intensity of the 
scattering in a given direction. Changing the ratio will change the profile shape, 
as shown in Figure A.7 (m/o = 10), and Figure A.8 (m/o = 20). The shadowing 
effect is neglected and only far-field is calculated. The Isakovich model calculates 
the dispersion of waves for random surface [Ref. 3]: 


т |1 — cos v cos E, cos Ó + sin Û sin E, 
с | (sin № + ѕіп E,)2 

1 уту2 cos*w + cos’? Е, – 2со5 соз В соз 0 
е. ESAE 2” 


p = 
а 


(4.8) 


а 
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Figure A.5: Bullington Model — Sine Square Profile. 


This model is much more complicated than the others, since it has two more 
parameters: 
Е, 15 the observation angle, and, 


0 1s the azimuth angle. 


In order to illustrate results comparable to those for the other modes, we 
assume that the azimuth angle is zero (@ = 0). For simplicity, we calculate the 
reflection coefficient for two cases: zero observation angle, E, — 0? (transmitter 
and receiver at the same level), and E, = —30? (transmitter higher than the 


receiver ). 
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Figure A.6: Davies Model. 
Using the assumption 0 = 0, we simplify Equation (A.8) to: 


m | — cos cos E, Lx Es] 


Pese T g (sin + sin E,)? 


| 1 sm\2 соз? + cos? FE; — 2 cos) cos E, 
е5 Е | ) (sin V + sin E,)? (4.9) 


а 


For the case of zero observation angle (E; = 0°), Equation (A.9) becomes even 


more simple: 


m / 1 — cos 1 /mNx2 cos? y + 1 — 2 cos 
ЕЕЕ ЕСІН” 


sin? y 8 Хо sin? у 


This case is illustrated in Figures A.7 and A.8 for two different values of 


т/с. For —30° observation angle Equation (A.11) applies. The result is shown in 


Figure A.9. 


р 


REFLECTION COEFFICIENT 


_ m 2 — ,/3cos ) — sin 
P&--9 7. | S(siny- 0.5)? 


E (ау соз?) + 0.75 — V3 cos % 
‘Exp = o a и ——ə— —nsy— > D. — —. 


Е | (4.11) 
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Figure A.7: Isakovich Model with Zero Observation Angle (£; = 0°), 


па о = 10, 
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Figure A.8: Isakovich Model with Zero Observation Angle, (Е; = 0°), 
mjo— 20. 
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Figure A.9: Isakovich Model with Observation Angle of —30°,m/o = 20. 
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APPENDIX B: COMPUTER PROGRAM LISTING 


This appendix gives the listing of the computer program described in Chapter 
VII. RUN is the management function, which calls the other functions in order. 
DATAIN is the function responsible for the input data for further calculations. 
STARTPRO is the main function which does all the major calculations, as well as 
creating the look-up tables, and comparing the ratios of the actual received signals 
with the calculated ratios. SIM is the simulation function, which simulates the 
received signals from the missile. RHFIND is the function which translates the 


location of the missile in the table into a range and height. 


xx RUN 
[0] КОМ 
C1] A MANAGMENT FUNCTION 
222 "ЕМТЕН ЅЕД 5ТАТЕ:' 
[3] 1 1 
CH] '0-CALM 1-SMOOTH 2SDIGH 
ГЕЯ 13-MODERATE 4 -ROUGH S-VERY ROUGH! 
[6] '6-HIGH 7-VERY HIGH! 
[7] 1 í 
[8] SEA<QO 
[9] a ENTERING ROUGHNESS PARAMETER 
[10 1 1 


[11] DATAIN 
[12] a INITIALIZATION 
[13] ЕР:Е 5ТАЕТРЕО НЕЗ 
[1uJ Rñ STARTING 
>[15] +Y(FLAG=1 )/STO 
[16] 05 END OF SESSION IF RANGE CANNOT BE FOUND 
[373 RHFIND 
[psu GO TO CHECK REAL RANGE AND DETCH 
> (1907 570550 


96 


xxDATAIN 


[0] DATAIN 

(13 а INITIALIZATION FUNCTION 

[2] а 00000000 

[3] MINHEIGT+4 

ОЯ а 'MINHEIGI' IS THE MINIMUM TX. HEIGHT (LESS 1 METER) 
55] A 00000000 

[6] HEICTST-36 

7 а 'HEIGTST' IS THE NUMBER OF HEIGHT STEPS 

[8] n 00000000 

(93 HEIGTRES<1 

10] A 'HEIGTRES! IS THE RESOLUTION OF HEIGHT STEPS 

ші) а 00000000 

ET? ) С<0.975 

EIS] a 00000000 

bit] р+1.025 

#15] а 'C!,'!D!' ARE LOWER AND UPPER LIMITS FOR OK (STARTPRO) 
mic] A 00000000 

5177 RMIN-1000 

EIS) a 'RMIN' IS MINIMUM RANGE 

[19] RMINS-RMIN 

201 а 'RMINS' IS MINIMUM RANGE FOR SIMULATION 

HOT] а 00000000 

[22] NSTEPS-125 

623) а 'NSTEPS' IS THE NUMBER OF RANGE STEPS 

[24] а 00000000 

[25] RESO<100 

[26] а 'RESO' IS THE RANGE STEPS RESOLUTION (LENGTH OF EACH STEP) 
027] о 00000000 

[28] HNSTEPS-130 

[29] а 'HNSTEPS' IS THE NUMBER OF HIGH RESOLUTION RANGE STEPS 
[301 а 00000000 

[31] HRESO+<50 

N22] A 'HRESO' IS THE RANGE STEPS (HIGH) RESOLUTION FOR SMALL 
[33] а RANGES 

[ uj а DDDDDDDD 

35] SIMSTEP<50 

[36] а IGIMSTEP! IS THE BANGE STEPS RESOLUTION IN THE SIMULATION 


NI7] a FUNCTION (VSIM) 
[38] F<(1999+(?17001))+1000 


[39J n 00000000 

[40] aF<9.35 

EET) a 'F! TS FREQUENCY (IN GHZ); 2718 GAZ 

[42] ERP-10000000x?100 

[43] a 00000000 

[uu ] aERP«-100000000 

(005 ] а 'ERP! IS TRANSMITTER EFFECTIVE RADIATED POWER (IN WATTS); 
[46] а 1E7+1E9 WATT (1004120 DBM) 

[47] RI«20000x(500*12300)21000 

(u8] a 'RI' IS INITIAL RANCE (IN METERS); 107 0ā00:168,000 METER 
(u9] а 00000000 
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E50) 
IE» 
[52] 
6531 
[Su] 
D 
[551 
S7 
(581 
(59) 
[60] 
(611 
(621 


[63] 
C64) 
[65] 
C66] 
C67] 
C68] 
[69] 
[793 
671) 
E72 
[73] 
Cu] 
-[75] 


ARI<15400 

HAT<U+736 

A 09000000 

ачаже? 

а 'HAT' IS TRANSMITTER HEICHT (IN METERS); это 
nRnHR3<1u+5716 

A 'HR3' ARE RECEIVER'S ANTENNAS HEICHTS (IN METERS) 

а 'HR3'! IS BETWEEN 15730 МЕ 

a 00000000 


HR3- 12 1" 16 18 20 

' RANGE. INTERVAL FOR SIMULATION: ! 
'MINIMUM RANGE (IN METERS): І, ФБЕМІК 
‘MAXIMUM RANGE (IN METERS): 


', SAMIN+(NSTEPSXRESO )+(HNSTEPSXHRESO ) 


'SEA STATE:! 

SEA 
A 'SEA! IS ROUGHNESS PARAMETER 
' FREQUENCY (IN GHZ):' 

E 

'INITIAL RANCE (IN KILOMETERS! 


RI*1000 : 

'"TRANSMITTER HEIGHT (TN МЕЛЕН) =ЛЕЛОЛЕ И 
НАТ 

(RECEIVER ANTENNAS HEIGHTS (IN METERS):! 
НЕЗ(АНЕЗ] 
A DISPLAING THE INITIAL SIMULATION DATA 
>0 
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**STARTPRO 
[O] ПИПИН ИЕ КОЛА НАТО РА ТИК ЛИВРА?ЕТІ;ЕТ?2ТЕТЗ?ЕТИЕТ5;ЕТ12:ЕТІ1З; 
ПИ а oc eon Soto Su: SITS: R12:R13:R14:R15:0K12: 

Ore ORM ТО ЛЕМРОТЕИРІ З ТЕМР2: 586 дрі 5р2: Врз :ярч: 

ОЕ ДИ DELTA? DELTAS: DELTAY:DELTAS:TETA1:TETA2: 
БЕЛДЕР ЛЕЛ 757 KII STO TKI AIL LIS?RABOI;REHO2:RAHO3 


ЕНОЧ;ЕНО5 
Г 1] з MAIN CALCULATION AND COMPUTATION FUNCTION 
[2] [oc 
23 1 Renin tanec OxXtaNoTEPS ), CHRESOXHNSTEPS )+RESOXINSTEPS 
СУ] a EU G ICE VECTOR (UHNSTEPS' STEPS OF ‘HRESO' METER 
c5] A EACH, ЛО МОТЕРБ О STEPS OF  !RESQU METER EACH), 
(61 д STARTING AT 'RMIN' 
[7] НА=МТМНЕТСТ+НЕТСТЕЕ5хХ\НЕТСТ5Т 
(81 A HA TS TRANSMITTER HEIGHT VECTOR ('HEICIST' STEPS OF 
[9] а 'HEIGTRES' METER EACH) 
#10) GRH<(HEICTST , (NSTEPS+HNSTEPS ) )p1 


Ël1] а IGRH' IS MATRIX WHICH SHOWS ALL THE LOCATIONS OF 
(121 а CANDIDATE RANGE AND HEIGHT COMBINATIONS 


[13] LAMBDA<0.3#F0 

C14) a 

15] Бвр1е(((НА-НЕ[1] )х2 )9. +6«%2 )#0.5 
[16] бр о (((НА-НЕ [2] )«2 )9. +Е«2 )#0.5 
817] RD3«-CCCGIA-HRL3] )*2)9. *R«2)*0.5 
[18] RDu«-CCCOHA-HRUu] )*2)9. *R«2)*0.5 


[19] RDS<( ((HA-HRES] )*2)°.+R*2)*0.5 
[#20] а 'RD' IS THE DIRECT DISTANCE BETWEEN TRANSMITTER 
21] а AND RECEIVING ANTENNAS с EL) 


[22] РЕГТА1+(2хНАхНЕ[1))°. +В 
[23] DELTA2<(2xHAXHR[2] )o. +R 
C24) РЕГТАЗ е(2 ХНАХНЕ [3] )о ,+Е 
[25] рЕрТтаЧе(2хХНАХНЕ [4] )9.:8 
[26] DELTA5«(2xHAxHR[5] )°.+R 


27] a КрЕБРГА!" IS THE DIFFERENCE BETWEEN THE DIRECT 
[28] а PATH AND THE REFLECTED PATH FOR EACH ЕХ. 
[29] A ANTENNA (RESPECTIVELY) 


[30] ТЕТД1 19 (9 (ФрЕРр1)р]НЕ(1]-н4)еЕрі 
[31] TETA2- 10(R(OpRD2)pl| HRUL2J -HA) * RD2 
[32] TETA3< l1o(&(@$pRD3)p|HR[3J-HA)+RD3 
[33] ТЕТАЧ 19 (9(ФоЕрч)ріне[ч]-на)зЕрч 
[34] TETA5«- 10(R(OpRD5)pl HRE5]-HA)*RD5 


[35 ] A 'TETA' IS THE OBSERVATION ANGLE BETWEEN THE RX. 
[36] д ANTENNA AND THE TRANSMITTER (FROM THE LATEST 
F 571 а POINT OF VIEW)- FOR EACH RX. ANT. (RESPECTIVELY) 


[383 XI1<(Q(@$pRD1)pHR[1J+HA)+RD1+DELTA1 
[39] XI2-(R(OpRD2)pHRU2]^*HA)*RD2*DELTA2 
[40] XI3-(R(ÓpRD3)pHRL3]*HA) * RD3*DELTA3 
[u1 ] XI4<(Q(OpPRDY )pHRO4]+HA )t+RDY+DELTAY 
uro ] XIS«-(R8(OoRD5)oHRUÜS1*HA)*RDS*DELTAS 
[43] A "ХІ! IS THE SINE OF THE GRAZING ANGLE FOR EACH 


(45) а RECEIVING ANTENNA (RESPECTIVELY) 
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>[45] 
>[u6] 
>[u 7] 
-[48] 
-(49] 
-[50] 
> [51] 
>52) 
(53] 
[54] 
(ЕБ Д 
[55 
[575 
C58] 
>[ 59] 
[60] 
[61] 
[62] 
+ [63] 
C64] 
C65] 
[66] 
>[67] 
[68] 
[69] 
ШІСІ 
»(711 
72] 
[73] 
[74] 
э [75] 
[76] 
[774 
[78] 
>[79] 
(80] 
(81] 
(82] 
-(83] 
C384] 
[85 
(86] 
(87] 
(88] 
(89] 
(90] 
[91] 
[92] 
[93] 
C94] 


+(SEA=0)/S0 
Э СОВА TJS 
XOSES-O p 52 
+(SEA=3)/S3 
>(SEA=4)/S4 
+(SEA=5)75S5 
> (БЕД-6)/55 
(СЕЛ ЕЛ УУЛ? 
a 
SO: SH-SHL«-SHH«O 
>СОМТ 
S1:SH+<0,032S5 
SAL<0 
5НН<0.065 
-CONT 
5.2 оке 375 
585-0.065 
SHH<0.21 
CONT 
53%254<0,265 
SAL<0.21 
5НН<0.32 
>СОМТ 
S4:SH<0.43 
SHL<0.32 
SHH<0.54 
-CONT 
O9:5H«*0.7 
SHL<0.54 
5НН<0.86 
-CONT 
S6:SH<1.08 
SHL<0.86 
SHH<1.3 
>СОМТ 
$7:5Н=1.75 
SULETI 
SHH<2.2 
-CONI 
a "SHL AND 


ARE THE LOWEST AND HIGHEST 


a (RESPECTIVELY) RMS WAVES HEICHT FOR EACH SEA STATE 


CONT: 


RHO1<* 2x((02)XxSHxXIA*LAMBDA)*x2 
ВНО2+ж 2х ( (02 )хбАхХІ2+ГАМВрА)х 2 
RHO3<* 2x((02)xSHxXI3+LAMBDA )x2 
РНОч+ж 2x((02)xSHxXIU* LAMBDA )x*2 
RHOS<* 2x((02)xSHxXIS+LAMBDA )x2 
n 'RHO' IS THE REFLECTION COEFFICIENT FOR EACH 


a RX. ANTENNA 
a AMENT MODEL 


(RESPECTIVELY) - CALCULATED AFTER 
(SEE VSIM) 
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[95] ET1<RHO1~=(10(( (02 )x ((AAx2 )+H2A01]xHA))°.+(LAMBDAxR))) 
C96] ET2<2HO2x(10(( (02 )x(CHA*2 )+HROT2]xHA))°.+(LAMBDAXxR))) 
(9871 ET3-RHO3x(10(((o02)xCCHAx2) *HRÜS]xHA))9.* CLAMBDAxR))) 
(98] ETu-RHOux(10(CCOo2)xCCHÀX*2) *HRUG])xHA))9. * CLAMBDAXxR))) 
poo ETS-RHOSx(10(((02)xCCHA*2) *HRUS1xHA))9. * CLAMBDAxR))) 
(100] a "ЕТ1!:!ЕТ5! ARE MULTIPATH ATTENUATION FACTORS FOR 
ШІСІ) A EACH RECEIVING ANTENNA, RESPECTIVELY (MATRIX) 
102] ЕТ12<ЕТ1:ЕТ2%1Е713 
(1031 ЕТ13<ЕТ1%:ЕТЗ%1Е713 
Сид ЕТ1Ч<ЕТ1+ЕТЧ+1Е 13 
(1051 ЕТ15<ЕТ1%:ЕТ5%1Е713 
(1061 а 'ET12'+'ET15' ARE RELATIONS BETWEEN THE MULTIPATH 
61021 -A FACTORS FOR EACH PAIR OF RX. ANTENNAS, RESPECTIVELY 
(1081 а (МАТВІХ) 
(1091 RCV:SI+«SIM 
МӘ) а 'SI! ARE RECEIVED SIGNAL LEVELS FOR ALL THE RECEIVING 
ШЕШІ) а ANTENNAS (VECTOR) 
12] REL: 
[113] 5112<(рСЕН9)р51111:51(21%1Е713 
[118] 5113<(рСЕН9)р51111%51ІГЗ31І%1Е 13 
B1i5] SI1u-(oCRH)oSIU[1]*SI([Uu]-*1E 13 
(116] SI15-«(oCRH)oSI(11*SIU(S]-^1E 13 
ЕНІП а (5112!:15115! Д2Е RELATIONS BETWEEN RECEIVED SIGNALS 
[118] a PAIRS (MATRICES) 
[119] COMP: 
[120] Е12<5112%ЕТ1241Е713 
1121] Е13<5113:ЕТ13%1Е713 
[122] R1u<SI1u+ET18u+1E 13 
[123] Е15<5115%ЕТ15%1Е 13 
(24 а 'R12'+'R15' ARE RATIO MATRICES BETWEEN ACTUAL 
шоў aA AND THEORETICAL VALUES 
[126] OK12-(R122C)^(R12sD) 
[127] OK13<(R132C)A(R13SD) 
F 128] OK1u<(R1u2C)A(R1usD) 
[129] OK15S<(R152C)A(R15SD) 
БІЗІ А 10К121%!0К15! ARE LOCATIONS OF RECEIVED SIGNALS 
(41311 A WHICH ARE EQUAL TO THE PRE-CALCULATED VALUES 
[132] a (ACTUAL RELATIONS EQUAL TO THE THEORETICAL ONES) 
[01033] aA (MATRICES) 
[134]  MATCOMP: 
[135] CRH«OK12^0K13^0K1*4^OK15 
(136] a 'CRH'! ARE LOCATIONS WHICH ARE EQUAL ALSO BETWEEN 
1137) а ANTENNAS PAIRS 
+>[138] -(1-*/,CRH)/RHFIN 
(139] a CHECK IF THERE IS A SINGLE COMBINATION OF RANGE 
[1450] а AND HEIGHT 
> (141) +(FLAG=1)/RHPIN 
[142] +RCV 
(143) a IF MORE THAN ONE COMBINATION, GO TO ANOTHER RECEIVING 
шц) а SESSION (AS LONG AS THE RANGE IS HIGHER THAN THE 
[155] а MINIMUM RANGE). IN CASE OF A SINGLE COMBINATION CO TO 
C146] a TRANSLATE INTO REAL READING (RANGE AND HEICHT) 
+(147] RHFIN:+0 
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ж«51М 


Co] A<+SIM+ LAMBDA +A1:A2:A3- AU: AS - SYP Shee see 
Ста а RECEIVED SIGNAL AHPLITUDES SIMULATION САШЕ 
(23 a RECEIVING ANTENNA), INCLUDING TRANSMITTER 
C3] а MOVEMENT INFLUENCE 
Cu] LAMBDA«-0.3*F 
C5] FLAC<0 
C6] SXI<«"30((|#R3-HAT)=RAI) 
tom а "SAIT IS THE CRAZINC ANCER 
-(8] +(SEA=0)/LAB1 
(9] SHR<SHL+(1? (SHH-SHL )x1000 )#1000 
Сао ж 'SHR' IS THE RMS WAVES HEIGHT FOR SIMULATION 


> Л >LAB2 
C198 LAB1:SHR<0 
(131 LAB2:SRHO<x 2x(( (02 )xSHRxSXI+LAMBDA ) )*2 


(181 в 

C15] A THE SIMULATION USES AMENT MODEL FORCAL CUTATINE 

C16] а THE REFLECTION COEFFICIENT (p): 

[21735 а 

(151, aA 

C19] а р =ЕХР. (=2х(2хРІХН xSIN. CXI )#LAMDA )x*2 ) 

СаО а S RMS 

21] +A 

C22] “A 

[23] 
A1<SRHOC1)x (ERPx4x (100 ¢ (02 )x (HR3 01) +HAT )xHAT )+(LAMBDAXRI )) )x*2 )+ 
o4 )xAI «2 

C24] 
A2-SRHOC(2J]xCERPxux(10(C(02)x (HR3U2] *HAT)xHAT) * (LAMBDAXxRI)))*2)*€ 
ou )xRIx2 

2c 
A3<SRHO[3]x (ERPxux(1o(((o2)x(ER3[3]+EAT)xRHAT)+ (LAMBDAxRI1 D sS 
ou )xRIx2 

(26] 
Au«-SRHO(uJxCERPxux(10(((02)x (CHAZ C4] tHAT XAT ) + (LAMBDAXxRI)))x2)*( 
o4 )xRI x2 

[27] 
А5<=#$ЕНО{5]хХ(ЕВРхцх(1о0(((о2)хХ(НВ3{51]1+НАТ)хНАТ)+(ЬАМБРАхЕТ)))х2)+( 
ou )xRI*2 


(28] a '"SRHO' IS THE ‘REFLECTION COEFFICIENT Bio K 

(291 а EACH RECEIVING ANTENNA (RESPECTIVELY) 

359 A+A1,A2,A3,A4,A5 

За а tAtIS А VECTOR OF RECEIVED AMPLITUDES (IN EACH 

(321 а ANTENNA-RESPECTIVELY) 

Us ЛЕ О Ка 

(351 а UPDATING THE TRANSMITTER RANGE (ІМ '‘SIMSTEP' METER Stee 
ра + (RIZRMINS)/FRET 

(36] t 1 

Ca] 'END OF SIMULATION SESSION?! 

(387 FLAG+1 

(391 a RANGE HAD NOT BEEN FOUND YET, BUT IS LESS THAN THE 


(u0] a MINIMUM RANCE, OR END OF SIMULATION SESSION 
-(u1] RET:+0 


x * RHE IND 
[0] RHFIND:INDH R n C МА ГАК :ТУБУЕС: БОС: УЕССЕН 
s] a FINDING AND DISPLAYING THE RANGE AND HEIGHT 
Вага) a OF IHE TRANSMITTER 


ЕЭ INDH«-100000xMINHEIGT-*HEIGTRESXxiHEIGITST 

СЧ) а ' INDH' IS AN HEIGHT INDEX 

[5] INDR-EMIN*(HRESOxiHNSTEPS),(HRESOxHNSTEPS)* (RESOxiNSTEPS) 
C6] а 'INDR' IS A RANGE INDEX 


[574] MATHR«-INDH9.*INDR 

{8 ] а 'MATHR' IS A RANCE / HEIGHT INDEX MATRIX 

[9] INDVEC+,MATHR 

C10] FIND: 

E11] VECCRH-,CRH 

ОДЕ CHANGING MATRICES INTO VECTORS 

E13] LOC-VECCRH/INDVEC 

[14] а SELECTINC THE CORRECT LOCATION OF THE TRANSMITTER 
Е(51 E OC! TS THE BX. LOCATION) 

C16] Е<100000|20С 

yaa) H-LLOC*100000 

[18] n TRANSLATING INTO REAL RANGE ('R!) AND HEIGHT ('H') 
выкаў 1 | 

E207] t! 

[ 24] 'TRANSMITTER RANGE (IN KILOMETERS):' 

#22) R+1000 

£23) ' TRANSMITTER HEIGHT (IN METERS):! 


[24] Н 
[25] ЈЕ. 
[26] ni 
[271 'UPDATED ACTUAL RANGE (IN KILOMETERS):! 
[28] (SIMSTEP+1000)+RI+1000 
>[ 29] +0 
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